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ABSTRACT 

As part of the SPLASH survey of the Andromeda galaxy (M31) and its neighbors, we have obtained 
Keck/DEIMOS spectra of the compact elHptical (cE) satelhte M32. This is the first resolved-star 
kinematical study of any cE galaxy. In contrast to previous studies that extended out to r < 30" ^ 
1 rf ~ 100 pc, we measure the rotation curve and velocity dispersion profile out to r ~ 250" and 
higher order Gauss-Hermite moments out to r '--^ 70". We achieve this by combining integrated-light 
spectroscopy at small radii (where crowding/blending are severe) with resolved stellar spectroscopy at 
larger radii, using spatial and kinematical information to statistically account for M31 contamination. 
The rotation curve and velocity dispersion profile extend well beyond the radius (r ^ 150") where the 
isophotes are distorted. Unlike NCC 205, another close dwarf companion of M31, M32's kinematics 
are regular and symmetric and do not show obvious sharp gradients across the region of isophotal 
elongation and twists. We interpret M32's kinematics using three-integral axisymmetric dynamical 
equilibrium models constructed using Schwarzschild's orbit superposition technique. Models with a 
constant mass-to-light ratio can fit the data remarkably well. However, since such a model requires 
an increasing tangential anisotropy with radius, invoking the presence of an extended dark halo may 
be more plausible. Such an extended dark halo is definitely required to bind a half-dozen fast-moving 
stars observed at the largest radii, but these stars may not be an equilibrium component of M32. The 
observed regularity of the stellar kinematics, as well as the possible detection of an extended dark 
halo, are unexpected if M31 tides are significant at large radii. While these findings by themselves do 
not rule out tidal models for cE formation, they suggest that tidal stripping may not be as significant 
for shaping cE galaxies as has often been argued. 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — galaxies: individual (M32, 
NCC 221) — galaxies: spectroscopy — galaxies: local group 



1. INTRODUCTION 

Andromeda's (M31) nearest companion, M32, is 
our closest example of a compact elliptical (cE), a 
rare dwarf galaxy type. Including M32, there are 
fewer than t en cE galaxies known within 100 Mpc 
dRool [19651: IKing fc Kised ITflTl ls"mith CasteUi et al. 
120091: iChilingarian et al.l l2009t iChilingarian fc Bergonc 
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120101 ). Compact elliptical galaxies have luminosities 
of ~ 10^ comparable to dwarf elliptical galaxies 
(dE), but with significantly smaller effective radii (ros ~ 
0.1-0.7 kpc), leading to notably higher surface bright- 
nesses (lieS ^ arcsec" IChoi et al.l [200l 

iMieske et all [20051 : IChilingarian et al.ll2009D . While cEs 
do appear to be a rare galaxy type, the discovery of 
these objects has been slow due to their elliptical galaxy- 
like appearance at distances beyond the Local G roup 
(jZiegler fc BendeH[T998HDrinkwater fc Gregei[l998[ ). 

Compact elliptical galaxies are consistently found 
in close projectio n to a massive parent g alaxy (5 < 
''proj < 80 kpc, IChilingarian et all l2009( ) . indicating 
that gravitational effects play some role in their evo- 
lution. One formation scenario proposes that cEs 
are remnant cores of tidally stripped 'normal' galax- 
ies ('Faber"1973'; 'Nieto"1990; 'Bckki et al."2001'; 'Graham| 
U002; Chilingarian & Bergond 2010; Huxor ct al. 2010!), 
although the rarity of cEs relative to normal elliptical 
and spiral galaxies suggests this does not happen of- 
ten, or that cEs are short-lived. Another formation the- 
ory suggests that cEs are low-mass classical elliptical 
galaxies that were either captured by or formed in the 
poten tial well o f a m assive neighbor (jNieto fc Prugniel 
Il987t iBurkcrt, 119941 ). This latter scenario is sup- 
ported by the position of cEs on the fundamental plane, 
at the low l uminosity end of the classical elliptical 
galaxy trend (jWirth fc Gallagheii[l98l lKormendviri985l: 
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iNieto fc Prugnie]||1987D . If these objects are indeed low 
mass normal elliptical galaxies, then their rarity implies 
a steep fa U off at the faint end of the galaxy lu minosity 
function (jBinggeli. Sandage fc TammannI 119881 ) . How- 
ever, recent observational improvements indicate that 
the light profiles of many cEs, including M32, are better 
fit by a two component r^/" bulge + exponential profile, 
characteristic of a disk galaxy, rather than the histori- 
cally used r^/^ law, typical of elliptical galaxies (Graham 
[2002). 

M32 is the nearest example of the cE class and an ex- 
cellent specimen for examining cE properties. As is the 
case for other cEs, M32 lies close to its parent galaxy at a 
projected separation of only 22' (5 kpc) from M31's cen- 
ter. Our relative proximity to M32 allows us to resolve 
individual stars in this object — at least for outer radii 
where stellar crowding is less severe. Photometric evi- 
dence suggests that M32's physi cal distance from M 31 is 
similar to its projected distance: iChoi et aD (120021 ) show 
that distortions in M32's outer elliptical isophotes are 
consistent with the hypothesis that these two galaxies 
are tidally interacting. Tidal interactions are a possible 
explanation for the galaxy's unusual s tellar population 



gradient and li ght profile (|F abcr 19 731: lO'Connelll 
Davideel 120001: (Pavidge et al.i i2000i :" lBekki et all 
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Grahaml 120021 ). Integrated light spectroscopic studies 
suggest the presence of a younger, mor e met al-rich stel- 
lar population at the center of M32 (iRo se et al. .2005i : 
iCoelho. Mendes de Oliveira fc Cid Fernan dcs 200!|, and 
references therein) which is perhaps the result of a tidally 
induced nuclear starburst (|Bekki et al.ll200il ). 

Previous studies of M32's internal kinematics have 
primarily focused on the inner regions. The steep 
rotational velocity gradient at the center of M32 
indicates the presence of a central dark mass. High 
resolution imaging and integrated light spectroscopy 
has confirmed the presence of a central black hole with 
Mrh ~ 2-4 X 10^ Goodman fc Led [1981 
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I2010D. While these studies 



have been appropriate for determining the mass of the 
central dark object, their limited radial extent (generally 
r < 10") does not provide much information on the 
wider galaxy environment. Past attempts to study the 
kinematics of M32 at lar ger radii (r < 60") have pro- 
duced conflicting res ults (lTonrvlll984t iCarter fc JenkinsI 
I1993D . Specifically, [Ibnol (|198ljf ound the velocity 
disper sion to increase outwards while ICarter fc Jenkini 
([1993') found it to be decrease outwards. Both authors 
profess problems in their measurements due to compli- 
cated sky subtraction and, in particular, contamination 
from M31. 

In this paper, we use a combination of integrated light 
and resolved stellar spectroscopy to obtain an accurate 
kinematical profile of M32 out to ~ SvcS — a much larger 
radius than has been previously possible. In the inner re- 
gions of M32, crowding is so significant that individual 
stars cannot be resolved, while in the outer regions, the 
integrated light spectrum becomes very noisy due to the 
presence of M31 and the steep light profile of M32. This 
is the first attempt to combine spectra of individual stars 



with integrated light spectroscopy to obtain a complete 
picture of M32. This work is part of the SPLASH Sur- 
vey (Spectroscopic and Photometric Landscape of An- 
dromeda's Stellar Halo), aimed at the study of M31 and 
its satellites. 

This paper is organized as follows. In §[2l we measure 
the internal kinematics of M32 out to a radius of ^ 8reff 
1 kpc) using resolved stellar light and ~ 3res (~ 
0.3 kpc) using integrated light. We then present velocity, 
velocity dispersion and Gauss-Hermite moment major- 
and minor- axis profiles from these measurements and 
compare to previous measurements. Using these results, 
we construct an axisymmetric, three- integral model in §[3| 
to obtain an estimate of M32's mass and M/L. Finally, 
in §[4|we summarize and discuss our results. 

Throughout this paper, we assume that M32 has the 
same distance modulus as M31: (m — M )n — 20.54 ±0.03 
(785 ±25 kpc. lMcConnachie et al.l[2005l ). At M32/M31's 
distance, 1" is equivalent to 3.8 pc. While the distance 
to M32 has been measured ind ependently (jJensen et al.l 
l2003t iKarachentsev et al]|2004f ) , it is consistent with the 
distance to M31. Moreover, the M32 distance estimate 
is affected by crowding problems and M31 contamina- 
tion. The lack of obvious signs of dust extinction in 
M32 suggests that the satellite lies in front of M31's disk 
(jFord. Jacobv fc JenneHll978f) . but a precise M31-M32 
distance is yet to be established. 

2. STELLAR KINEMATICS OF M32 

In this section, we give a detailed account of our M32 
kinematical measurements. In ii l2.1l we present measure- 
ments from the resolved stellar kinematics. In ii l2.2l we 
present measurements from integrated light. For read- 
ers not interested in these details, an integrated view of 
M32's kinematics and a comparison to previous studies 
are presented in ii l2.3l 

2.1. Resolved Stellar Population Spectroscopy 

In this subsection, we discuss our observational setup, 
including photometric and astrometric measurements 
(i) l2.1.ip . mask placement (ii l2.1.2p . identification of iso- 
lated sources (i) l2.1.3p . target selection and slit mask de- 
sign (i) l2.1.4p . and observing details (i) l2.1.5l) . We then 
discuss the reduction and analysis of the spectra, includ- 
ing the data reduction process (ii l2.1.6p . measurement of 
individual stellar velocities (§ l2.1.7p . quality assessment 
(i) l2.1.8p . and detection of and velocity measurements 
for serendipitously detected stars (§[27L9|). FinalUy, in 
ii l2.1.10l we make velocity and velocity dispersion mea- 
surements along M32's major and minor axes via maxi- 
mum likelihood analysis. 

2.1.1. Pre-Imaging 

We derive photometric and astrometric catalogs of 
stars in M32's general vicinity from archival data ob- 
tained with the MegaCanF^ imager on the 3.6 m Canada- 
France-Hawaii Telescope (CFHT). These data are in the 

http : //www . cf ht .hawaii . edu/Instruments/Imaging/MegaPrimeT] 

Megafrime/iVlegaCam is a joint project ot (Jt'HT and 
CEA/DAPNIA at the Canada-France-Hawaii Telescope (CFHT) 
wliicii is operated by the National Research Council (NRC) of 
Canada, the Institut National des Science de I'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, 
and the University of Hawaii. The observations were obtained by 
the MegaCam instrument team in November 2004. 



M32 Stellar Kinematics 



3 



form of 2° x 2° mosaic images in the g' , r' , and i' bands 
centered on M31. 

The mosaic image in each band was built from four 
separate MegaCam pointings each with a field of view of 
1° X 1°. Each MegaCam pointing consists of five dithered, 
slightly overlapping exposures with integration times of 
5 X 45 s in each of the g' and r' bands and 5 x 90 s in the 
i' band. Thus, the effective integration times are 225 s in 
each of the g' and r' bands and 450 s in the i' band, r' 
bands, and 450 s in the i' band. 

The data were obtained under photometric conditions. 
Individual exposures were detrended using CFHT's data 
reduction pipeline. Elixir, and stacked into mosaic im- 
ages using the SWarp software (jBertin et al\ f2002]. The 
pixel scale on the mosaic image is 0'.'36, a factor of 2 
coarser than the native MegaCam pixel scale. The me- 
dian seeing FWHM is 1'.'2, I'.'O and 0'.'7 in g', r', and i', 
respectively. In order to highlight point-like sources, a 
high-pass filtered image is constructed by subtracting a 
smoothed version of the mosaic image from the original. 

For the purposes of this M32 study, we analyze a 36' x 
36' section of the i' and r' mosaic images centered on M32 
(given the high degree of crowding in this region, the g' 
image is not particularly useful due to its relatively poor 
seeing and the r' image is only used in the construction 
of Figure |31 neither the g' nor r' images are used for 
spectroscopic selection) . The size of this image section is 
chosen to comfortably allow for the placement of multiple 
Keck/DEIMOS multislit spectroscopic masks (each mask 
covers about 16' x 4'), as shown in Figure [TJ 

We generate star lists over the high-pass filtered r' 
and i' MegaCam mosaic images using t he FIND task 
in the DAOPHOT photometry package ()Stetsonl [1991 
and carry out preliminary aperture photometry of all 
detected sources. A spatially varying point spread func- 
tion (PSF) template is iteratively derived from a set of 
bright, isolated stars whose neighbors have been sub- 
tracted using these star lists. This PSF template is then 
fit to all sources in the catalog for each frame using the 
ALLSTAR module to produce accurate photometric cat- 
alogs. Best fit PSF templates are then subtracted from 
the high-pass filtered image to create a subtracted image 
that shows the residuals due to imperfect PSF subtrac- 
tion and missed sources. The source find procedure and 
photometry procedures are then repeated on the PSF 
subtracted images several times to identify faint objects 
missed in earlier passes, and this yields a final photo- 
metric catalog of over 10*^ stars. Next, the magnitudes 
are roughly calibrated based on the tip of the red giant 
branch (TRGB) in each filter; precise photometric cali- 
bration is not needed for this project since the purpose of 
our photometric catalog is solely to assign rough i' mag- 
nitude ranges from which to select and set priorities for 
spectroscopic targets. As a final step, we astrometrically 
correct all of our x and y positions by computing trans- 
formations based on the USNO A2.0 catalog. The final 
RA and DEC of all sources are accurate to ~0.2". We 
demonstrate our photometric characterization of sources 
in the i' image by illustrating both original and PSF sub- 
tracted images in the vicinity of M32 in Figure [2] 

FigureOshows the i' and r' photometry for stars within 
the 5% membership probability contour centered on M32 
(see Figure [1] right). Note that the color- magnitude 
diagrams (CMDs) in Figure [3] are limited by the poor 



quality of the r' data, and are therefore used here for 
illustration purposes only. The plots on the left show 
the photometry of targets in the velocity range of M32 
(-275 <v < -125kms"\see i) l2.1.10|) . while the figures 
on the right plot the remaining non-M32 like targets. All 
four CMDs look very similar, thus illustrating the diffi- 
culty of photometrically preselecting stars that are likely 
to be M32 members. At best, the CMDs can be used 
to preselect stars at or below the TRGB (i' > 20), and 
above our spectroscopic limit («' < 22). 

2.1.2. Mask Placement 

The two largest obstacles to spectroscopy of the re- 
solved stellar population of M32: contamination by M31 
and crowding. We consider each of these in turn. 

M32 is projected against the high surface brightness 
inner regions of M31. The probability of a star being an 
M32 member can be estimated by the M32/M31 surface 
brightness ratio at that location. The two-dimensional 
surface brightness of M32 is mod eled as a s er ies o f 
concentric ellipses based on the iChoi et all ()2002f ) 
measurements of /-band surface brightness, ellipticity 
and position angle as a function of radius. The two- 
dimensional surface brightness of M31 is modeled as the 
sum of a Sersic bulge and an exponential disk based 
on published V-band su rface brightness data and a 
mean colo r V - I = 1.6 CWaherbos fc Kennicutt! 119871: 
prit chet fc van den Bergh 1994; Guhathakurta et al.l 
120051 llrwin et al.l 1200511 . The 5% and 50% 
M32/(M31-hM32) /-band surface brightness ratios 
are shown in Figure [T] (for a more detailed mapping see 
Figure [5]) . The probability of M32 membership drops 
off very rapidly with increasing radius because of M32's 
steep surface brightness profile. 

Near the center of M32, however, it is difficult, if not 
impossible, to find isolated spectroscopic targets in the 
high-pass filtered CFHT image. Even the most luminous 
RGB stars (that are otherwise ideal for spectrosco py) ar e 
badly blended. Based on visual inspection (see ij |2.1.3p . 
we avoid the inner rM32 ^ 1' region, corresponding to a 
surface brightness of fii = 18.9 mag arcsec^^. 

The arrangement of our five Keck/DEIMOS multislit 
masks is shown in Figure [TJ Each mask covers « 16' x 4'. 
The masks provide coverage of the rM32 < 150" region 
where M32's isophotes are regular, and the r M32 > 150" 
regio n where isophotal distortion is observed (jChoi et al.l 
I2002| ). The first slitmask, M32_l, is centered on M32 
with the long axis rotated to a position angle of 160°, 
i.e. the approx imate posit ion angle of the inner ellip- 
tical isophotes (|Choi et al.l [20021. The remaining four 
multislit masks, M32_2— M32_5, are oriented to optimize 
coverage of the outer regions. M32's steep brightness pro- 
file implies that there is only a limited region in which 
neither contamination nor crowding is too severe; our 
arrangement of masks ensures good coverage of this re- 
gion, while also covering M31's inner spheroid and disk 
for other SPLASH survey science. 

2.1.3. Identifying Isolated Sources 

Our chosen placement of the five Keck/DEIMOS 
masks defines the footprint over which spectroscopic tar- 
gets are selected. We next identify stars within this 
footprint which are least affected by crowding/blending. 
We use two criteria to reject targets: (1) cases where 
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Figure 1. Positions of the five observed Keck/DEIMOS multislit masks (M32_l - M32_5, black rectangles) shown within a 36' X 36' 
section centered on M32. Each DEIMOS mask is ~ 16' X 4' Left: Masks overlaid on the Choi et al. (2002) KPNO Burrell Schmidt B-band 
image. Red circles mark the location of our spectroscopic targets. The dashed white circular contour at the center of M32 denotes the 
i' surface brightness magnitude limit at which crowding becomes so significant that individual targets can no longer be resolved in the 
CFHT/MegaCam image (/^j/ ~ 19 mag/arcsec^, rM32 ~ !')• The surface brightness gradient seen in this image (brighter towards the 
north-west) is contaminating light from M31. Right: Masks overlaid on the high-pass filtered CFHT/MegaCam image used for target 
selection. The dashed black line marks the location where M32's isophotes begin to show distortion l(Choi et al.ii2002r i. The elliptical-like 
blue and red contours mark the location of M32's 5% and 50% membership probabilities, respectively, determined on the basis of model 
fits to the 2D surface brightness distributions of M31's bulge/disk and M32. 

DAOPHOT finds one or more bright neighbors tliat are 
close to, but distinct from, the target (we refer to this as 
crowding), and (2) cases where an apparent single source 
in the DAOPHOT catalog is a poor fit to the PSF (we 
refer to this as blending). 

We address the issue of crowding by rejecting a target 
star if it has a neighbor in the DAOPHOT catalog that is 
so close/bright that the PSF of the neighbor significantly 
overlaps that of the target. Based on visual inspection 
of crowded regions of the image, we have come up with 
an empirical criterion. Any target that has even a single 
neighbor satisfying the following relation is ehminated 
from the list of potential spectroscopic targets: 



itgt 



nbr 



Ftgt,nbr| 

0.8" 



- 3.0 



(1) 



where /tgt and /nbr are the apparent magnitudes of the 
target and neighbor, respectively, and rtgt,nbr is the pro- 
jected distance between the two objects. Of the sources 
in the DAOPHOT catalog in the magnitude range / = 
20 — 22 (t he range used to select spectroscopic targets, 
see ? l2.1.4p . « 10% pass this crowding test; the surviving 
fraction increases with target brightness over this mag- 
nitude range. 

We address the issue of blending by visually inspecting 
the images at the locations of the stars that survive the 
crowding test. This inspection includes both the high- 
pass filtered and PSF-subtracted versions of the I'-band 
CFHT/MegaCam image, as illustrated in Figure[2l Each 
target is fiagged as unblended (high priority), marginally 
blended (medium priority), or badly blended (low prior- 
ity), depending on the degree to which its image resem- 
bles the PSF on the high-pass filtered image and the 
strength of systematic residuals at its location on the 
PSF-subtracted image. 



These de-blending exercises are only good to a point 
as we are limited by the O'.'S seeing (FWHM) of the 
CFHT/MegaCam image. The seeing FWHM was signif- 
icantly better than this during the Keck/DEIMOS spec- 
troscopic observations (§ l2.1.5p . As a result, further de- 
blending of sources is carried out in the sp atial and spec- 
tral domains, as discussed in ii l2.1.9l See iDormanet al.l 
(|201^ ] for a discussion of an automated procedure for 
identifying blended sources in this data set. 

2.1.4. Mask Design 

Targets are prioritized for spectroscopic observation 
based on two criteria: level of blending and magnitude. 
The first prioritization, based on level of blending, di- 
vides the targets into three lists: list 1 - unblended, list 
2 - marginally blended, and list 3 - badly blended (as 
discussed in ii l2.1.3p . The second prioritization, based 
on magnitude, assigns priorities wi thin each list based 
on M 31's TRGB magnitude of 20.5 (jMcConnachie et all 
120051 ): the highest priority targets (priority 1) have I = 
20.5 — 21, and the lowest priority targets (priority 2) have 
/ = 20 - 20.5 or J = 21 - 21.5. Targets with magnitudes 
outside the range / = 20 — 21.5 are excluded. The distri- 
bution of targets across these magnitude ranges is shown 
in Figure [H 

Five Keck/DEIMOS multislit masks are designed us- 
ing A. C. Phillip's dsimulator softwareF^. The software 
takes as input the multiple target lists (lists 1-3) orga- 
nized by target priority (priorities 1-2). Each w 16' x 4' 
Keck/DEIMOS mask is populated with targets from list 
1, in order of decreasing priority, followed by list 2, and 
so on. Our selecti on process is identical t o tha t discussed 
in Appendix A of iGuhathakurta et al.l (|2006f l. with the 

^■^ [http : //www ■ ucollck ■ org/ ~phllllps/delmos_ref /masks ■ html | 
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Figure 2. Example of the CFHT/MegaCam image quality in various regions. Left: Outline of the M32_l DEIMOS mask (grey) on the 
high-pass filtered CFHT/MegaCam image with three 0'.'7 X I'.'O blue boxes defining representative areas zoomed into on the right. Middle: 
The column of plots demonstrates the quality of the image and the density of targets at the various locations on the mask. The selections 
show: (top) region near M32, and (bottom) the southern part of the mask. The red circles mark the location of the sources detected in 
the photometry. Right: The PSF subtracted images. 



following modifications: target prioritization, minimum 
distance between target and slit end (1'.'65), and dis- 
tance between adjacent slitlets (CK'S). The location of 
the slitlets selected for observation is shown in Figure [TJ 
The five masks contain a total of 883 slitlets. 

2.1.5. Observations 

Five multislit masks in the region of M32 were ob- 
served between November 20 07 and August 20 08 using 
the DEIMOS spectrograph (|Faber et alj I2003D on the 
Keck II 10 m telescope. The arrangement of the masks is 
shown in Figure[T] The observing details are summarized 
in Table m 

All five multislit masks were observed with the 1200 
line mm~^ grating. This configuration yields a spa- 
tial scale of 0'.'12 pixel" "'^ and a spectral dispersion of 
0.33 A pixel"^. We set the central wavelength to 7800 
A, corresponding to wavelength range of ~ 6450 — 9150 
A. The exact wavelength range for each slit varies as a re- 
sult of location on the multislit mask and/or truncation 
due to vignetting. The wavelength region is chosen to 
target several spectral features including the strong Ca 
II triplet absorption feature present in RGB stars. The 
anamorphic distortion factor for this grating and cen- 
tral wavelength is 0.606. Therefore, the 0'.'8 wide slitlets 
subtend 4.1 pixels. Better still, excellent seeing condi- 
tions (~ 0'.'6) during observations provide somewhat bet- 



ter spectral resolution yielding an average resolution of 
3.1 pixels = 1.0 A FWHM. 

Useful spectra are obtained from 786 of the 883 slitlets 
(89%). The success rate would have been even higher 
were it not for the fact that w 25% of the slitlets on mask 
M32_3 were lost due to buckling of the mask during its 
insertion into the DEIMOS focal plane at the time of 
observations; fortunately the buckling took place at the 
ENE end of the mask away from M32. 

2.1.6. Data Reduction 

The five Keck/DEIMOS multislit masks are processed 
using the spec2d and spec Id software (version 1.1.4) de- 
veloped by the DEEP Galaxy Redshift Survey team at 
the University of California, BcrkelejO- Briefly, the re- 
duction pipeline rectifles, flat-fleld and fringe corrects, 
wavelength calibrates, sky subtracts, and cosmic ray 
cleans the two-dimensional spectra, and extracts the one- 
dimensional spectra. We give a more detailed description 
of the reduction process below. 

First, the reduction pipeline rectifles curved spectra 
into rectangular arrays by applying small shifts and in- 
terpolating in the spatial direction. One-dimensional slit 
function, two-dimensional flat-field and fringing correc- 
tions are then applied to the spectra. The wavelength 

|http : //astro ■ berkeley . edu/~ cooper/deep/spec2d/ 1 
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Figure 3. Color-magnitude diagram (CMD) based on the i' and r' photometry from the high-pass filtered CFHT MegaCam image (for 
illustration only due to the poor quality of the r' data). Stars shown arc within the 5% M32 membership probability contour (see FigurefTI 
right) and outside 1' of M32's severely crowded center (see Figure [T] left). The absolute calibration of the data is approximate and 
based on the TRGB magnitudes for M31 from McConnachic ct al. (2005). The black circles/squares denote Kcck/DEIMOS spectroscopic 
targets with measured velocities. The numbered rectangular outlines mark the location of our spectroscopic priorities: (1) primary and 
(2) secondary. Top Left: Black circles show M32 candidate members (i.e. stars with velocities in the range —275 < v < —125km s~^) 
within 1' < rM32 < 2.'5. Top Right: Black squares show M32 non-members (i.e. stars with velocities in the range v < —275 km s~^ or 
V > —125 km s~^) within 1' < rM32 < 2f5. Bottom Left: Black circles show M32 candidate members with riv[32 > 2.'5 (but within M32's 
5% probability contour, i.e. rM32 < 5f3). Bottom Right: Black squares show M32 non-members with rM32 > 2.'5 (but within M32's 5% 
probability contour, i.e. rM32 < 5.'3). Note the similarity between all four populations, thereby illustrating the difficulty of photometrically 
disentangling these two populations. 



solution of the rectified spectra is obtained by fitting a 
polynomial to the arc lamp lines (precise at the O.OlA 
level). The two-dimensional spectra are then sky sub- 
tracted and cosmic ray cleaned. Sky is identified by col- 
lapsing the two-dimensional spectra in the wavelength 
direction to locate spatial positions along the slit that 
are least affected /unaffe cted by targets and serendipi- 
tous sources (see ^ 12.1.91 for a discussion of serendipitous 
sources). Each two-dimensional spectrum is sky sub- 
tracted by fitting a B-spline model (wavelength as a func- 
tion of two-dimensional position: x, y) to the night sky 
emission lines in the baseline portion of the spatial inten- 



sity profile. This careful sky subtraction is of particular 
importance around the Ca ii triplet region due to the 
presence of bright night sky lines; poor sky subtraction 
would reduce our ability to accurately measure stellar 
velocities. The two-dimensional exposures are then com- 
bined along with cosmic ray rejection and inverse vari- 
ance weighting to create a single mean two-dimensional 
spectrum for each slit. 

Last, the target is identified and its one-dimensional 
spectrum extracted. Targets are located on the two- 
dimensional spectrum by identifying the peak bright- 
ness distribution in spatial intensity profile, obtained by 
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Table 1 

Keck/DEIMOS Multislit Mask Exposures 



Mask 
Name 


Observation 
Date 


a 

{J2000.0) 


5 

{J2000.0) 


P.A. 

(deg) 


tcxp 

(m) 


Seeing 
FWHM 


No. of 
Slits 


No. of Useable 
Target Velocities 


No. of Useable 
Serendip Velocities 


M32_l=' 


2007 Nov 14 


00 42 38.3 


+40 51 34.0 


160 


2 X 20 


0'.'5 


194 


189 (97%) 


73 


M32_2 


2008 Aug 03 


00 43 03.8 


+40 55 07.7 


70 


3 X 20 


o'.'e 


184 


169 (92%) 


17 


M32_3 


2008 Aug 03 


00 43 11.6 


+40 52 34.7 


-110 


3 X 20 


0'.'7 


191 


134 (70%)*^ 


10 


M32_4 


2008 Aug 04 


00 42 13.9 


+40 54 44.2 


105 


3 X 20 


o'.'e 


143 


137 (96%) 


117 


M32_5 


2008 Aug 04 


00 42 13.9 


+40 52 02.6 


-75 


3 X 20 


0'.'6 


171 


157 (92%) 


81 


Total: 














883 


786 (89%) 


298 



Note. — Units of right ascension (a) arc in hours, minutes and seconds. Units of declination {5) are in degrees, arcminutes and 
seconds. 

^ The "M32_l" mask was originally named "M32" at the time of submission of the mask design. 

^ Buckling of the M32_3 DEIMOS multislit mask at time of observation adversely affected sr; 25% of its slitlets. 

bv iSimon fc Gehal (|2007[ ). The scale factor 1.85 is de- 
termined from duplicate radial velocity measurements of 



collapsing the two-dimensional spectrum in the wave- 
length direction. Target one-dimensional spectra are ex- 
tracted from the two-dimensional spectra using a small 
spatial extraction window centered on the target. The 
one-dimensional spectra are re-binned into logarithmic 
wavelength bins with 13.8km s~^/pixel. The final result 
is a wavelength calibrated, sky subtracted, cosmic ray 
cleaned one-dimensional spectrum for each target. An 
illustration of this process is shown in Figure |3^. 

2.1.7. Cross Correlation Analysis 

Line-of-sight (LOS) velocities for resolved targets are 
measured from the one-dimensional spectra using a 
ISimon fc Gehal (l200l modified version of the visual in- 
spection software zspec, developed by D. Madgwick for 
the DEEP Galaxy Redshift Survey at the University of 
California, Berkeley. The software determines the best- 
fit LOS velocity for a target by cross correlating its 
one-dimensional science spectrum with high S/N stel- 
lar templates in pixel space and locating the best-fit in 
reduced space. The ten best-fit templates, LOS ve- 
locities, reduced-x^ values and cross-correlation errors 
are reported. The stellar templates used in the cross- 
correlation analysis cover a wide range of stellar types, 
F8 to M8 giants, subgiants and dwarf stars, and metal- 
licities, [Fe/H] = -2.12 to +0.11 (Simon & Gcha 2007). 
The observing setup for the templates is nearly identical 
to that discussed in § 12.1.51 with the exception that tem- 
plate stars are observed with 0'.'7 wide slitlets, the min- 
imum slit length is set to 4" (to allow for adequate sky 
subtraction), and the template stars are trailed across 
the slit. 

A-band telluric corrections and heliocentric corrections 
are calculated and applied to the measured LOS veloci- 
ties. The A-band telluric corrections, which account for 
velocity errors associated with the slight mis-centering 
of a star in a slit, arc determined using the method dis- 
cussed in lSohn et a l. (2007). 

LOS velocity errors are determined for each star by 
scaling the cross-correlation based velocity error us- 
ing duplicate radial velocity measurements of stars ?, 
K. Gilbert, private communication The average LOS ve- 
locity error for M32 RGB The velocity error for each star 
Av is estimated to be: 



Av = V(1.85 X At;cc)2 +2.22 



(2) 



stars. The average LOS velocity error for M32 RGB stars 
is 4km s~^. 



Each 



2.1.8. Quality Assessment 
two-dimensional spectrum, one-dimensional 



where Avcc is the cross-correlation based error and 
2.2 km s~^ is the systematic velocity error as determined 



spectrum, and corresponding Doppler shifted template 
match, are visually inspected in zspec and assigned a 
quality code Q based on the reliability of the fit. This 
process allows the user to judge the quality of a spec- 
trum and reject instrumental failures and poor quality 
spectra. Velocity measurements based on two or more 
strong spectral features are assigned (5 = 4 (82% of tar- 
gets). Velocity measurements based on one strong fea- 
ture plus additional marginal features are assigned Q = 3 
(6% of targets). Spectra that contain no strong features, 
low S/N and/or instrumental failures are assigned Q — 2 
(8% of targets). For cases in which zspec did not return 
an accurate velocity measurement, but visual inspection 
of the one-dimensional spectrum showed an obvious ve- 
locity shift, the velocity is manually marked and assigned 
Q = 1 (< 1% of targets). Foreground stars used for the 
purpose of alignment are assigned Q — —1 (3% of tar- 
gets). Additional details on this quality code assignment 
can be found in Guhathakurta ct al. ( 2006() . 

2.1.9. Serendipitous Sources 

Upon visual inspection of the one-dimensional and two- 
dimensional spectra during the quality assessment phase 
outlined in § 12.1.81 some fraction of the slits clearly show 
that the full length of the slitlet intersects more than 
one star: the target star and one or more serendipitously 
detected stars, known hereafter as a serendips. These 
detections occur frequently in our target region due to 
the severe crowding and blending in the CFHT Mega- 
Cam data. It is common for the full length of the slit 
to intersect multiple isolated sources in addition to the 
target (i.e. neighbors) as a result of the severe crowd- 
ing in the region. Sources that are nearly blended with 
the target are also commonly found; this is in part be- 
cause the de-blending exercises for target selection (see 
S I2T31) are good only to the 0'.'8 seeing limit of the CFHT 
MegaCam data. The better angular resolution of the 
Keck/DEIMOS spectroscopic data ('^ 0'.'6) allows for the 
spatial resolution of stars that can not be resolved in the 
CFHT MegaCam data. 

Serendips are detected via one of two methods: 
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(a) Target Only 

Ha 
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(c) Composite Target — Consists of Two Superposed Sources 
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Figure 4. Examples of target stars and serendipitous sources (serendips). Each panel (a-c) shows a slitlet's two-dimensional spectrum 
(top left) where the wavelength axis x runs horizontally (zoomed in to highlight specific spectral features), and the spatial axis y (position 
along the slitlet) runs vertically. The horizontal streak(s) within each two-dimensional spectrum is the stellar continuum from RGB 
candidate(s). The spatial location of continuum is determined by collapsing the two-dimensional spectrum in the wavelength direction 
and identifying peaks in the intensity profile (top right, black). One-dimensional stellar spectra (bottom, black) are extracted from the 
two-dimensional spectrum using a small spatial extraction window (top, red) that straddles the continuum of the desired star. The LOS 
velocity of a star is measured by cross-correlating Doppler shifted stellar templates with the one-dimensional stellar spectrum until a best 
match is found (bottom, blue). Note that the best- fit Doppler shifted stellar template (blue) is offset from the stellar spectrum (black) for 
the purpose of illustration, (a) Example of a slitlet intersecting a single star (i.e. the target star), (b) Example of a slitlet intersecting 
two spatially resolved stars (i.e. the target star and an off-target serendip). The one-dimensional spectrum and best-fit stellar template 
shown correspond to the serendip. (c) Example of a slitlet intersecting two spatially blended stars (i.e. the target star and an on-target 
serendip). The presence of multiple stars is revealed by the two sets of absorption features in the one-dimensional spectrum. 



through continuum detections that are offset from the 
primary target in the spatial direction (referred to as 
off-target serendips) , or by the detection of spectral fea- 
tures that are offset from the primary target in the spec- 
tral direction (referred to as on-target serendips). In a 
couple rare cases, we detect stars using a combination 
of the two methods (referred to as off-target superim- 
posed serendips). We discuss the details of these detec- 
tion methods in turn. 

Off-target serendips are visually identified as addi- 
tional brightness peaks in the spatial intensity profile 
that are offset from the target and spatially coincident 



with spectral continuum in the two-dimensional spec- 
trum (see Figure [JJd). This includes cases where the 
spectral continuum is distinct from the target and in- 
stances where it is partially blended with the target (i.e. 
the spatial intensity profile shows a peak with an asym- 
metric wing). Once the location of the off-target serendip 
has been visually identified in the spatial intensity pro- 
file, an extraction window is manually placed on the two- 
dimensional spectrum and the one-dimensional spectrum 
is extract ed. T he reduction process then proceeds as out- 
lined in ii l2.1.7l - ? I2.1.8I the zspec software is run to find 
the best template match to the one-dimensional spec- 
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trum, the LOS velocity is measured, telluric and helio- 
centric corrections are applied, the LOS velocity error 
is calculated, and a quality is assigned. LOS velocities 
are measured for 244 off-target serendips over the five 
Keck/DEIMOS multislit masks. 

On-target serendips are cases where the spatial inten- 
sity profile shows only a single peak (i.e. looks like one 
star) but where 2 sets of spectral absorption features 
are evident (Ca ii, TiO, etc., see Figure |4]:). The min- 
imum velocity separation for which we detect two dis- 
tinct superimposed velocities is Av ^ 50 km s~^ (which 
is slightly greater than the ~ 35 km s^^ FWHM of our 
instrumental resolution near the Ca ii triplet). We define 
the fainter of the two superimposed stars as the on-target 
serendip. The best template and LOS velocity match to 
the on-target serendip is determined during the quality 
assessment phase. In cases where the ten best-fit solu- 
tions reported by zspec include a mix of fits for both 
the target and the on-target serendip, the template and 
LOS velocity best matching the on-target serendip is se- 
lected from the list and assigned a quality code {Q =3, 
4) . In cases where the solutions reported by zspec do not 
show any good matches to the on-target serendip, the 
best-fit is determined by manually shifting stellar tem- 
plates until a good match is found; these fits are assigned 
(5 = 1. Once the LOS velocity has been determined, a 
heliocentric and an average telluric correction is applied, 
and the LOS velocity error is calculated. LOS veloci- 
ties are measured for 52 on-target serendips over the five 
Keck/DEIMOS muhislit masks. 

Off-target superimposed serendips are a combination of 
the two categories discussed above. They consist of cases 
where continuum that is offset from the primary target 
turns out to be two perfectly blended serendips at the 
same spatial location. The one-dimensional spectrum of 
the off-target superimposed serendip is extracted via the 
method outlined for off-target serendips. The zspec soft- 
ware is then run on the one-dimensional spectrum, and 
analysis proceeds as for on-target serendips. We find 2 
such instances of off-target superimposed serendips over 
the five Keck/DEIMOS muhislit masks. 

LOS velocities are measured for a total number of 298 
serendips in all three classes. The number of serendip 
LOS velocities measured per mask is summarized in Ta- 
ble ffl 

2.1.10. Maximum Likelihood Analysis of the Line-of-Sight 

Velocities 

A subsample of the measured target and serendip 
LOS velocities is selected for further kinematical analysis 
based on their probability of M32 membership, as deter- 
mined by sky position. Stars with a probability P > 5% 
(see Figure [1]) are selected as potential M32 candidates. 
This reduces the stellar sample from 1,084 to 482 stars 
and corresponds to a radial range of 0.2 — 1.4 kpc. In or- 
der to determine if strong velocity and/or velocity disper- 
sion gradients are present along M32's major and minor 
axes, the subsample of stars is further divided into eight 
subregions: four quadrants divided into two probability 
ranges. The four quadrants, outlined in Figure[5l are cen- 
tered on M32's axes as follows: north-north-west (NNW) 
major-axis (—65° < 4> < 25°, spanning 0.2 — 0.9 kpc), 
south-south-east (SSE) major-axis (115° < < 205°, 
spanning 0.3 — 1.4 kpc), west-south-west (WSW) minor- 



axis (205° < < 295°, spanning 0.2 - l.Okpc), and east- 
north-east (ENE) minor-axis (25° < < 115°, span- 
ning 0.3 — 1.0 kpc). The lack of symmetry between the 
contours in M32's SSE and NNW quadrants is due to 
differing amounts of M31 light contamination; the M31 
contamination is significantly worse in M32's NNW re- 
gion resulting in contours that extend further out on the 
SSE side of the galaxy. This effect can also be seen along 
the minor-axis of the galaxy where the contours extend 
further out on the ENE side of the galaxy. The two prob- 
ability ranges are defined so that each subregion contains 
a reasonable number of M32 stars: an inner region de- 
fined by P > 50% (0.2 < a < 0.6 kpc), and an outer 
region defined by 5% < P < 50% (0.5 < a < 1.4 kpc). 
These eight subregions are illustrated in Figure [5] 

We perform maximum likelihood fits of Gaussians to 
the LOS velocity distribution of stars in the M32 region 
(all stars with P > 5%) and in each of the eight sub- 
regions. While the true shape of the velocity structural 
components in these regions may differ from pure Gaus- 
sians, the use of such models seems appropriate given 
their ability to characterize the mean velocity and veloc- 
ity dispersion, the small number of velocity points being 
assessed, and the absence of any definite physical model. 
The individual stellar velocity errors are not included in 
the maximum likelihood analyses, since these errors (^^4 
km s~^) are much smaller than the velocity dispersion 
of any structural component in any of these regions, and 
therefore contribute insignificantly to the maximum like- 
lihood errors. The Gaussian fits to the M32-like resolved 
stellar LOS velocities in each region and corresponding 
68% and 90% confidence limit errors are summarized in 
Table [2] and discussed below in turn. 

Figure IHl shows maximum likelihood fits of sums of 
Gaussians to the LOS velocity distribution of stars for 
the entire M32 region (aU stars with P > 5%). The ob- 
served LOS velocity distribution is well fit by the sum 
of three Gaussians (red curve). The narrow Gaussian 
centered at v — —196.9 ± 3.0km s"'^ with a width of 
a = 29.9 ± 2.9km s^^ (solid black curve) represents 
35.2^35% of the stars in the region and is consistent 
with the systemic velocity of M32 (u^f = -200 km s"\ 
iFalco et aD 119991 ). The two additional populations seen 
in this region are well fit by a broad Gaussian centered 
at V = —350.9 km s~^ with a width of cr = 153.0 
km s~^ representing 41.9% of the stars, and a nar- 
row Gaussian centered at v = —386.9 km with 
a width of cr = 35.0 km s~^ representing 22.9% of 
the stars; these two additional populations are consis- 
tent with the mean v elocity of M31's inn er spheroid 
u^^f 300km s-i (IGilbert et al.l 120071) and M3rs 

disk v^^^^ ^ —400 km s~^, respectively ([Dormanet al.l 
[20T1 . 

Figure [7] shows maximum likelihood fits of sums of 
Gaussians to the LOS velocity distribution of stars 
to each subregion located along M32's major-axis and 
minor-axis, respectively. In each subregion, the two 
Gaussians that represent M31's inner spheroid and disk 
(determined in the global fit to the M32 region) are 
held fixed while the best fit Gaussian parameters to M32 
members are searched. The mean velocity v, velocity dis- 
persion (T, and fraction iVM32/-^tot of M32 stars and the 



10 Howley et al. 



Table 2 

Maximum Likelihood Gaussian Best Fit Parameters to M32-like Resolved Stellar Velocities 



f RMS(r) ^ 




Av d 


Act f iVM32/A'tot S 


AAfM32/Aftot 1^ 


Confidence Limits: 




(68%) (90%) 


(68%) (90%) 


(68%) (90%) 



M32, All Quadrants, Inner + Outer 

174'.' 2 

Major Axis 



M32, NNW Quadrant, Inner 

-95"4 ll'.'g -198. 

M32, NNW Quadrant, Outer 

-201('2 43'.' 1 -2i 

M32, SSE Quadrant, Inner 

123'.'8 23'.'4 -n 

M32, SSE Quadrant, Outer 

225'.'6 

Minor Axis 



225'.'6 53^'5 -191.0 ±11-9 36.2 tfi tl^s° 0.311 



174'.'2 65'.'0 -196.9 ±3.0 +11 29.9 til -I'l 0-352 ±0.028 ILo^^H 



Qt;"4 11"Q 1QS7 +6-4 +10.8 r, . r, +6.0 +11.8 n 7A9 +0.092 +0.140 

-95.4 il. 9 -198.1' _g 5 _;^i 4 24.2 ^ q 0.^^42 -O.no -0.189 

9ni"9 ' zl'^"1 _99R zl +11-0 i +9.9 +19.3 n 1 Qfi +0.055 +0.094 

zoi. z "io. 1 _ii.4 -20.0 ■-'^•^ -7.4 -11.8 o.l»u -0.052 -0.084 

123'.'8 ' 23'.'4 ' -186.3 ±9.0 tlH 31.2 j:^;^ t^^^^ q ggg +0.107 +o.i65 



+0.080 



M32, ENE Quadrant, Inner 

inR"9 ISA 4 +7'0 +11.7 r,-, r, +7.8 +13.5 n 7HS +0.094 +0.139 

-106. 2 5. 1 -184.4 _g 4 3 26.2 2 -7.6 ^.788 -0.113 -0.195 

M32, ENE Quadrant, Outer 

-142('6 30'.'2 -204.6 t^jj 26.0 t^^g" 0.401 ±0.080 tf^m 

M32, WSW Quadrant, Inner 

119'.'3 21'.'4 -194.6 ±4.2 18.1 0.724 tg-O^I tg-.l^S 

M32, WSW Quadrant, Outer 

le^K"-? 9H"Q _17Q9 +11-7 +20-2 oq n +9-7 +18-7 907 +0.072 +0.120 

lOO. o zo. 3 _ii 3 -23.0 '-'^•^ -7.8 -21.8 -0.070 -0.113 



Note. — Results of the maximum likelihood Gaussian fits to the LOS velocities of the M32-like resolved stellar 
data (includes targets and serendips) in various quadrants located within M32's 5% predicted fractional light 
contribution contour (see Figure fsl. 

^ Median projected semi-major/minor axis distance from M32's center for the M32-like stellar population (i.e. stars 
with velocities in the range vicr) based on lChoi et al.l (1200211 /-band photometry. The full M32 region (all quadrants) 
and major axis quadrants (NNW and SSE) list the semi-major axis distance a. Minor axis quadrants (ENE and WSW) 
list the semi-minor axis distance qa. 

^ Root-mean-square of the semi-major/minor axis distance. 
^ Best-fit heliocentric LOS velocity. 

^ Error in best-fit heliocentric LOS velocity (68% and 90% confidence limits). 

Best-fit velocity dispersion. 
^ Error in best-fit velocity dispersion (68% and 90% confidence limits). 
^ Best-fit fraction of M32 stars based on kinematics. 

^ Error in best-fit fraction of M32 stars in the sample (68% and 90% confidence limits). 



corresponding 68% and 90% confidence limits for eaclr 
subregion are listed in Table [2l 

We compare the M32/M31 fraction predicted by our 
2D surface brightness models for the two galaxies to the 
fraction measured from the fit to the LOS velocity distri- 
bution in each subregion. While the two sets of fractions 
follow similar trends across the different subregions, the 
surface brightness based predictions of the M32 fraction 
tend to be systematically lower than the fractions mea- 
sured from the LOS velocity distribution analysis (but 
by only ^ 1 cr on average). In any case, this slight dis- 
crepancy does not affect any of the results of this paper. 
We only use the surface brightness based fraction predic- 
tions to define the boundaries of the subregions in which 
to carry out the kinematical fits, and these boundaries 
are fairly arbitrarily defined in any case. Nevertheless, 
we consider some possible explanations for this discrep- 



ancy between predicted and measured fractions. First, 
M31's disk surface brightness is not uniform across the 
M32 region and may well depart from our idealized 2D 
surface brightness model. Second, our spectroscopic tar- 
get selection tends to bias the kinematical sample against 
high surface brightness/crowded patches in M31's disk, 
resulting in a bias towards higher M32 fractions. Third, 
the translation from integrated T^-band surface bright- 
ness to RGB star count surface density is likely not the 
same for M32, M31's inner spheroid, and M31's disk, and 
may well vary with radius within the two galaxies. Fi- 
nally, there are uncertainties in our 2D surface brightness 
models for the two galaxies associated with /- to l/-band 
conversion of the surface brightness measurements. 

Our kinematically based measurement of the M31 con- 
tamination fraction is relevant for M3 2 stellar popula- 
tion studies. For example, in the recent iMonachesi et al.l 
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Inner: 0.25 <a< 0.55 kpc □ _ 
Outer: 0.48 <a< 1.38 kpc □ 
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Figure 5. A subsample of targets (closed circle) and serendips (open circle) used for kinematical analysis. Concentric open circles denote 
multiple serendipitous detections at a given location. The contours shown outline M32's predicted fractional /-band light contribution 
relative to M31's inner spheroid and disk. For the purpose of kinematical analysis, the M32 region is divided into four quadrants: NNW 
major-axis, SSE major-axis, WSW minor-axis, and ENE minor-axis. In order to determine if strong velocity gradients are present, each 
quadrant is further divided into three subregions based on the probability, P, of M32 membership (derived from two-dimensional surface 
brightness profiles of M31 and M32). The subregions for each quadrant are defined as follows: inner = 50% < P < 90% (corresponding 
to a semi-major axis distance of 0.25 ^ a < 0.55 kpc), and outer = 5% < P < 50% (0.48 ^ a < 1.38 kpc). These subregions are chosen 
(somewhat arbitrarily) so that the inner and outer subregions contain a reasonable number of M32 stars. Maximum likelihood fits of 
Gaussians to the LOS stellar velocities in each of these subregions are shown in Figure |6] 



12 



Howley et al. 




-600 -400 -200 

Velocity (km/s) 

Figure 6. Velocity histogram and maximum likelihood Gaussian 
fits to the LOS velocities of resolved stars located within M32's 
5% predicted fraction light contribution contour (see Figure [Sjl. 
The solid black curve shows the fit that is consistent with M32's 
systemic velocity = -200km s'^ [Falco et al.lll999l '). The 

dashed black curves have mean velocities and velocity dispersions 
consistent with a cold rotating component (M31's disk), with a 
hot non-rotating component (MSl's inner spheroid). The sum of 
the Gaussian fits is shown as a solid orange curve. The fits to the 
M31-like stellar population are held fixed during the analysis of the 
quadrants surrounding M32. 



()2011| ) analysis of M32's star-formation history from deep 
HST ACS/hRC CMDs, contamination by M31 stars is 
statistically accounted for using a control field whose lo- 
cation was chosen on the basis of the shape/orientation of 
MSl's disk isophotes. Our kinematical analysis provides 
an independent measurement of the M31 contamination 
fraction in the region of their study, albeit averaged over 
a larger area than the narrow HST ACS/HRC field. 

2.2. Integrated Light Spectroscopy 

We begin this subsection by describing observations, 
and end with measurements of the line-of-sight veloc- 
ity distribution (LOSVD) along M32's major and minor 
axes from the integrated light data. This subsection is 
outlined as follows. In ii l2 2.1[ we provide the observ- 
ing details. In ii l2.2.21 we summarize the data reduction 
process. In ij l2.2.31 '^6 make velocity, velocity dispersion 
and higher order Gauss-Hermite moment measurements 
along M32's major and minor axes. 

2.2.1. Observations 

Eight longslit exposures centered on M32 were ob- 
tained with Keck/DEIMOS between November 2007- 
2008. Six of the longslits were aligned with M32's major- 
axis {(j) = 160°), and two were aligned with M32's minor- 
axis (0 — 70°). Each longslit mask is ~ 16' long and de- 
signed with a series of 4 slits separated by small bridges 
in order to insure the structural integrity of the mask 
design. The observing setup for the longslit exposures 
is identical to that used for the multislit mask observa- 
tions discussed in ii l2.1.51 except that slit widths of O'.'S 
and I'.'O were used. The average seeing for these observa- 
tions was I'.'O, yielding an average spectral resolution of 
3.9 pixels — 1.3 A. The observing details for the longslit 
observations are summarized in Table [31 

2.2.2. Data Reduction 



Processing of the eight Keck/DEIMOS longslit masks 
from rectification through wavelength calibration is 
identical to the reduction procedure used for the 
Keck/DEIMOS multislit masks: the two-dimensional 
spectra are rectified, fiat-field corrected, fringe corrected, 
and wavelength calibrated using the spec2d reduction 
pipeline (see ? l2.1.6p . 

One-dimensional spectra are extracted at increasing 
spatial intervals from the center of M32 along the two- 
dimensional spectrum. The location of M32's center is 
determined by fitting a Moffat profile + 1st order poly- 
nomial to the intensity profile of M32 (obtained by col- 
lapsing the two-dimensional spectrum in the wavelength 
direction); the fractional spatial pixel location that cor- 
responds to the peak of the intensity profile fit is defined 
as M32's center. This information is used to convert the 
pixel positions along the two-dimensional spectrum into 
a distance from M32. One-dimensional spectra are then 
extracted using boxcar extraction windows with widths 
ranging from 1" to 5" (the actual size of the extraction 
window is determined by the signal-to-noise S/N of the 
region being extracted), and the Poisson errors are cal- 
culated. The one-dimensional spectra are re-binned log- 
arithmically in the wavelength direction into bins with 
13.8 km s^ /pixel. 

Next, sky subtraction is performed on the one- 
dimensional spectra. Extreme care is taken to prop- 
erly subtract the light from all contaminating sources, 
which includes M31, atmospheric air glow, and imper- 
fect fiat-fielding (resulting from differences in the illu- 
mination between the internal flat-field exposure and 
the on sky science exposure). Sky subtraction is per- 
formed separately for each wavelength bin. First, the 
intensity profile as a function of distance from M32's 
center is obtained for each wavelength bin. A normal- 
ized de Vaucouleurs' profile with — 29" is fit to and 
subtracted from each intensity profile outside the inner 
5"of M32's center, whe re excessive flux is present (jKenti 
Il987t iChoi et al.ll2002[ ): this produces a "M32-free" in- 
tensity proflle for each wavelength bin. Next, a 2nd or- 
der polynomial representing all the contaminating light 
sources is fit to each "M32-free" intensity profile at po- 
sitions beyond 90" from M32's center. The 2nd order 
polynomial fits are subtracted from the original inten- 
sity profiles, and then reassembled into one-dimensional 
sky-subtracted spectra. An example of the fitting and 
subtraction process is shown in Figure [S] 

Once sky subtraction is complete the major-axis spec- 
tra observed with the O'.'S slit width (see Table [3]) are 
Gaussian smoothed in the spectral direction in order to 
match the spectral resolution of the I'.'O slit width spec- 
tra (since the minor-axis science spectra were all observed 
with the same slit width they do not require smoothing 
for co-addition). The smoothing length (Tsmooth needed 
for the 0'.'8 major-axis spectra is determined to be 0.42 A 
using the following formula: 



Csmooth — m 



AA 

Ax 



!ln2 



(3) 



where W2 — I'-'O is the desired slit width, wi — 0'.'8 is 
the observed slit width, m = 0.606 is the anamorphic 
demagnification factor at the central wavelength 7800 A, 
AA = 0.32 A/pixel is the spectral dispersion, and Ax = 



M32 Stellar Kinematics 



13 



Velocity (km/s) 
-600 -400 -200 



Velocity (km/s) 
-600 -400 -200 



g L (a) NNW Inner 




. - (e) ENE Inner 




n 



10 

8 

6 

4 

2 


10 

8 

6 

4 

2 


10 

8 

6 

4 

2 


10 

8 
6 
4 
2 



-600 -400 -200 
Velocity (km/s) 



-600 -400 -200 
Velocity (km/s) 



Figure 7. Velocity histograms and maximum likelihood Gaussian fits to the LOS velocities of resolved stars along the major and minor 
axes of M32. Solid black curves show fits that are consistent with M32's systemic velocity [v^^-^ = —200km s~^. IFalco et al.1ll999 l). Dashed 
black curves (M31 components) have fixed ratios, velocities and velocity dispersions that are based on a global fit to the M32 region (Figure 
(BJl. The sum of the Gaussian fits is shown as a solid orange curve. Panels a-h show fits to the subregions outlined in Figure [5] (a) the 
inner region along M32's NNW semi-major axis, (b) the outer region along M32's NNW semi-major axis, (c) like (a) but for SSE major 
axis, (d) like (b) but for SSE major axis, (e) like (a) but for ENE semi-minor axis, (f ) like (b) but for ENE semi-minor axis, (g) like (a) 
but for the WSW minor axis, (h) like (b) but for the WSW minor axis. 

O'.'1185/pixel is the spatial scale. 

Next, the one-dimensional sky-subtracted spectra from 
individual exposures and original fine spatial bins are 
coadded, with cosmic ray rejection, into broader spatial 
bins. These broader bins are designed to achieve a mini- 
mum spectral S/N ratio per pixel of 25 averaged over the 
Ca II triplet region. This S/N threshold limits the radial 
extent of our kinematical analysis to < 90" and < 60" 
on the major and minor axes, respectively. About 15% 
of the individual spectra are excluded from the coadds as 



they were deemed, upon visual inspection, to suffer from 
bad subtraction of night sky emission lines and/or other 
systematic errors (e.g., slit edges and bad columns). In- 
clusion of these bad spectra in the coadds would lead to 
choppier, less symmetric kinematical profiles. No coad- 
dition is performed on spectra located at |rM32| < 2" as 
the individual spectra exceed the S/N threshold in these 
bright inner regions. 
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Table 3 

Keck/DEIMOS Long-Slit Exposures 



Mask 


Observation 
Date 


a 

( 19000 0^ 


s 

(J2000.0) 


P.A. 

(deg) 


tcxp 

(m) 


Seeing 
FWHM 


Slit 
Width 


Major.l 


2007 Nov 14 


00 42 41.87 


+40 51 57.2 


160 


3 


0'.'5 


1^'G 


Major_2 


2008 Nov 14 


00 42 41.87 


+40 51 57.2 


160 


3 


0'.'5 


l^'O 


Major_3 


2008 Oct 01 


00 42 41.87 


+40 51 57.2 


160 


5 


~ I'.'O 


I'.'O 


Major_4 


2008 Oct 01 


00 42 41.87 


+40 51 57.2 


160 


5 


~ I'.'O 


I'.'O 


Major_5 


2008 Nov 24 


00 42 41.87 


+40 51 57.2 


-20 


5 


I'.'S 


O'.'S 


Major_6 


2008 Nov 24 


00 42 41.87 


+40 51 57.2 


-20 


5 


1'.'3 


O'.'S 


Minor.l 


2008 Nov 24 


00 42 41.87 


+40 51 57.2 


70 


5 


1'.'3 


O'.'S 


Minor_2 


2008 Nov 24 


00 42 41.87 


+40 51 57.2 


70 


5 


1'.'3 


0'.'8 


Note. — 


The units of ri 


y\\t ascension {a 


) are in hours, minutes and seconds. The units of 



declination ((5) arc in degrees, arcminutcs and seconds. 



^ lO' 



B 10' 
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Figure 8. An example of sky subtraction and the night sky vari- 
ations in M32's integrated light profile. The plot shows light inten- 
sity (black diamonds) versus distance from M32's center in two 
wavelength bins. The red dashed line provides a fit to M32's 
light by normalizing a de Vaucoulours' profile, with r'f^ = 29" 
HChoi et al.ll200^ . to the observed intensities. The contaminat- 
ing light from M31, atmospheric air glow and corrections to the 
DEIMOS spectrograph response function is measured by subtract- 
ing off the normalized de Vaucouleurs' profile from the observed 
intensities and fitting a 2nd order polynomial to the residual inten- 
sities beyond 90" from M32's center (dashed blue line). The sum 
of the two profiles is shown as a solid black line. Top: Wavelength 
bin corresponding to night sky continuum. The gradient seen in 
the continuum (dashed blue line) is a result of light contamination 
from M31. Bottom: Wavelength bin corresponding to a peak of 
a night sky emission line. The M31 gradient is not visible against 
the dominant night sky emission. 



2.2.3. Measurement of Velocity, Velocity Dispersion and 
Higher Order Gauss- Hermite Moments 

The mean velocity v, velocity dispersion a, and Gauss- 
Hermite moments /i3 through hQ that best fit the ob- 
served LOSVD are determined using the p ixfit soft- 
ware developed and described by Ivan der Marel . (,1994 ). 



Briefly, the software determines the best-fit absorp- 
tion line strength parameter (7), LOS velocity (d), and 
velocity dispersion (cr) for each one-dimensional inte- 
grated light science spectrum by minimizing be- 
tween scaled, Doppler-shifted, Gaussian-broadened spec- 
tral templates and the science spectra (see Figure IHl- 
Symmetric and anti-symmetric deviations from Gaus- 
sianity are then measured by expanding the fits to in- 
clu de higher order Gauss-Hermit e moments (/i3 through 
fa6. lvan der Marel fc Fran"^ll993f ) for spectra with an av- 
erage S/N per pixel > 40. 

The rms uncertainties in our best-fit parameters x — 
[f , (7, /i4, ft,5, hQ] are determined using the Poisson 
error in the flux of the input science spectrum. Since 
the Poisson-based error estimates A-^x are bound to be 
underestimates and do not account for systematic er- 
rors such as imperfect subtraction of night sky emis- 
sion lines, template mismatch, and residual detector ar- 
tifacts, we attempt to derive more realistic error esti- 
mates, Aa; — fx A^x, where the error scale factor / 
is derived empirically by assuming that the true kine- 
matical profiles in the higher S/N outer regions of M32 
(|r| > 5", S/N > 40) are symmetric and smooth. The 
symmetry assumption implies that one half of the anti- 
symmetric profiles {v, h3 and hb) can be point-reflected 
onto its other half, and that one half of the symmetric 
profiles (cr, hA and h6) can be mirror-reflected onto its 
other half. The smoothness assumption implies that the 
folded data points a::data(''i) can be compared to a second- 
order polynomial Xpoiy(?') that has been fit to the data. 
Specifically, we require that deviations of the folded pro- 
files from the smooth polynomial has a per degree of 
freedom (xbof) '^^ unity: 



Xdof 



-E 



['^data (^i ) '^poly (^"i )] 



[Ax(r,)]2 

thereby defining the error scale factor: 

1 ■^-^ [a^data(''i) ~ ^Epoly (?'i)]^ 



1, 



{APx{n)\^ 



(4) 



(5) 



Error scale factors are calculated for each the major- 
and minor-axis velocity, velocity dispersion and com- 
bined Gauss-Hermite moment profiles. 

The LOSVD is analyzed over the wavelength range 
8470-8700 A (i.e., the region around the Ca 11 triplet: 
8498, 8542, 8662 A). This wavelength range is selected to 
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Figure 9. An illustration of the fitting of a composite stellar spectral template (red) to a 1" (4 pc) wide integrated light one-dimensional 
spectrum of M32. The composite template is constructed by optimally combining 16 stellar templates using a genetic algorithm. The three 
panels illustrate the steps used to measure the LOSVD of the M32 spectrum . (a) Observed M32 spectrum and an un-broadened composite 
stellar template (with no Doppler shift applied), (b) The Doppler velocity of the M32 spectrum is measured by shifting the un-broadened 
composite stellar template to match the features of the M32 spectrum, (c) The velocity dispersion of the M32 spectrum is determined by 
broadening the composite stellar template so that it matches the observed M32 spectrum. 

maximize the LOSVD "signal" while minimizing system- 
atic errors. A comparative analysis of the LOSVD dis- 
tribution between independent sections of the spectrum 
("blue": 6500-8470 A and "red": 8470-8900 A) indicates 
that systematic errors are relatively large in the "blue" 
portion. Scaling the Poisson -based errors on the LOSVD 
parameters to ensure "blue" vs. "red" agreement overpre- 
dicts the errors by at least a factor of two relative to those 
derived from the full spectrum and a scaling based on the 
assumption of smoothness and symmetry /antisymmetry 
in the radial LOSVD profiles. 

The spectral template used to analyze the LOSVD 
is constructed by combining weighted stellar templates. 
The choice of spectral template is of particular impor- 
tance as mismatch between the science spectra and spec- 
tral templa te can result in significant systematic errors 
(jRix fc Wh ite 1992). Given that M32 is composed of a 
variety of stellar types, it is not surprising that the galaxy 
spectra are not well fit by any single stellar spectrum. For 
this reason, a composite stellar template is constructed 
using a weighted linear superposition of stellar templates 
(details on the stellar templates can be found in ij |2.1.7|) . 
Since the stellar templates were observed with 0'.'7 wide 
slits, the templates are smoothed in the spectral direc- 
tion so that they match the resolution of the science spec- 
tra; the templates are smoothed by CTsmooth = 0.50 A for 
comparison to the I'.'O resolution major-axis data, and 
by CTsinooth = 0.27 A for comparison to the O'.'S reso- 
lution minor-axis data (see Equation [S]) . The optimal 
weights for the stellar templates are determined using a 



genetic algorithm (GA) as configured by iHowlev et al.l 
(,2008) . The GA locates the global minimum in the 
weight parameter space by minimizing the between a 
high S /N science spectrum and co-added weighted stellar 
templates, which have been Doppler shifted to match the 
science spectrum using cross-correlation and Gaussian 
smoothed using an initial velocity dispersion estimate. 
We assume that the velocity dispersions across the differ- 
ent stellar components are the same (i.e. not a function of 
age, stellar mass or metallicity) so that the co-added stel- 
lar templates are all smoothed with the same amount of 
velocity dispersion. Because template matching becomes 
increasingly difficult with lower S/N spectra, the proce- 
dure is run on a single high S/N spectrum (rM32 = — 2'.'0, 
S/N — 134) with a well measured velocity dispersion 
(o- = 56.2km s"\ Ivan der Marel et al.l[l994bl ). The in- 
tention is to use an unshifted, un-broadened version of 
the weighted composite stellar template to measure the 
velocity and velocity dispersion profiles for each science 
spectrum, thereby treating the abundance and popula- 
tion gradients in M32 as roughly constant with radius 
such that only the LOSVD, signal and noise of each spec- 
trum varies. Once the optimal combination of stellar 
templates is found, the continuum of the composite stel- 
lar template is adjusted to better match the science spec- 
trum. This is achieved by fitting a 5th order polynomial 
to the ratio of (science spectrum) / (broadened composite 
stellar template), both to the blue and red sides of the 
spectrum separately, and multiplying the resulting poly- 
nomial by the un-broadened template to remove any low 
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order frequency differences. Tlie final matcli between the 
science spectrum and composite stellar template is shown 
in Figures mjc-d) . 

Table |4] lists the best-fit parameters and scaled errors 
to the LOSVD as a function of radius. The error scale 
factors applied to the major-axis Poisson-based error es- 
timates are /„ = 3.8, fa = 2.8, and fh = 1.9. The error 
scale factors applied to the minor-axis Poisson-based er- 
ror estimates are fv = 2.1, = 1.9, and fh = 2.1. 
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Table 4 

M32 Integrated Light Profiles for Velocity, Velocity Dispersion, and 
Gauss- Hermite Moments 



r" dr^ v'= Ai;d a" Aa^ hz^ Aha^ Ha^ Ahi^ hs^ Ahs^ he^ Ahe^ S/N' Mask Namei 



M32 Major Axis (4> = 160°), Integrated Light 



-85^'0 


10"0 


-23.6 


7.3 


40.0 


7.9 


- 


- 


- 


- 


- 


- 


- 


- 


28 


- 


-TO'.'O 


5'.'0 


-19.5 


6.0 


40.6 


6.9 


- 


- 


- 


- 


- 


- 


- 


- 


33 


- 


-eo'.'o 


5'.'0 


-17.6 


5.0 


47.9 


5.1 


0.033 


0.085 


0.102 


0.102 


-0.002 


0.092 


-0.103 


0.101 


47 


- 


-50'.'0 


5'.'0 


-19.2 


2.5 


41.1 


2.8 


0.047 


0.058 


-0.067 


0.073 


-0.028 


0.062 


0.081 


0.071 


80 


- 


-40'.' 


5'.'0 


-18.9 


2.4 


42.6 


2.6 


-0.001 


0.052 


-0.040 


0.066 


0.053 


0.057 


0.038 


0.064 


81 


- 


-32'.'5 


2<'5 


-18.1 


2.1 


48.7 


2.4 


0.053 


0.037 


-0.019 


0.045 


-0.019 


0.041 


0.029 


0.045 


97 


- 


-irih 


2^'5 


-23.6 


1.3 


47.7 


1.4 


0.055 


0.024 


0.001 


0.029 


-0.005 


0.026 


0.017 


0.029 


151 


- 


-22'.'b 


2^'5 


-24.2 


1.4 


52.5 


1.5 


0.054 


0.021 


0.011 


0.026 


-0.002 


0.024 


-0.002 


0.025 


149 


- 


-iT.'b 


2^'5 


-30.1 


0.9 


51.8 


1.0 


0.048 


0.014 


0.008 


0.017 


0.002 


0.016 


0.010 


0.017 


222 


- 


~12'.'b 


2!.'b 


-36.2 


0.8 


54.4 


0.9 


0.048 


0.012 


0.017 


0.014 


0.022 


0.013 


0.001 


0.014 


250 


- 


-7'.'5 


2!.'b 


-43.8 


0.6 


55.7 


0.6 


0.050 


0.008 


0.016 


0.009 


0.019 


0.009 


0.002 


0.009 


360 


- 


-4'.'5 




-44.9 


0.9 


58.8 


1.0 


0.060 


0.011 


0.030 


0.013 


0.000 


0.012 


-0.015 


0.013 


240 


- 


-3'.'5 


0^'5 


-44.9 


0.7 


59.7 


0.8 


0.059 


0.009 


0.014 


0.011 


0.003 


0.010 


0.001 


0.011 


287 


- 


-2'.'5 


0^'5 


-45.7 


0.8 


60.8 


0.9 


0.058 


0.009 


0.024 


0.011 


0.006 


0.011 


-0.009 


0.012 


258 


- 


-I'.'g 


0('6 


-43.4 


1.1 


66.2 


1.2 


0.066 


0.011 


0.039 


0.014 


-0.009 


0.013 


-0.026 


0.014 


199 


Major_6 


-I'.'g 


o'.'e 


-44.9 


1.1 


62.8 


1.2 


0.067 


0.012 


0.034 


0.014 


-0.009 


0.014 


-0.019 


0.014 


199 


Major_5 


-1'.'3 


o'.'e 


-45.0 


0.9 


62.0 


1.0 


0.064 


0.010 


0.030 


0.013 


-0.006 


0.012 


-0.015 


0.013 


228 


Major_3 


-I'.'S 


O'.'e 


-45.3 


0.9 


62.2 


1.0 


0.067 


0.010 


0.028 


0.012 


-0.004 


0.012 


-0.010 


0.012 


236 


Major_4 


-o^'g 


o'.'e 


-42.8 


1.0 


70.1 


1.0 


0.069 


0.009 


0.014 


0.011 


0.002 


0.011 


-0.003 


0.011 


236 


Major_2 


-0'.'7 


o'.'e 


-22.6 


0.9 


81.3 


0.9 


0.056 


0.006 


0.010 


0.007 


-0.027 


0.007 


-0.009 


0.007 


312 
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-0"7 


O'.'e 


-23.7 


0.8 


79.1 


0.8 


0.062 


0.006 


0.015 


0.007 


-0.024 


0.007 


-0.011 


0.007 


319 


Major_5 


-OCl 


O'.'e 


-7.5 


0.7 


87.9 


0.7 


0.040 


0.005 


-0.017 


0.005 


-0.016 


0.005 


0.014 


0.005 


399 
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-oyi 


O'.'e 


-5.1 


0.7 


89.4 


0.7 


0.034 


0.004 


-0.016 


0.005 


-0.018 


0.005 


0.014 


0.005 


415 


Major_4 


O'.'S 


O'.'e 


14.7 


0.9 


87.3 


0.9 


-0.054 


0.006 


-0.012 


0.006 


0.029 


0.007 


0.009 


0.007 


322 


Major_2 


o^'s 


0'.'5 


22.9 


0.8 


81.6 


0.9 


-0.047 


0.006 


0.003 


0.007 


0.022 


0.007 


-0.004 


0.007 


320 


Major_6 


O'.'S 


0'.'5 


23.8 


0.8 


79.3 


0.8 


-0.056 


0.006 


-0.001 


0.007 


0.027 


0.007 


0.002 


0.007 


326 


Major_5 


I'.'G 


O'.'e 


44.3 


0.8 


66.2 


0.8 


-0.097 


0.008 


0.022 


0.010 


0.048 


0.009 


-0.012 


0.010 
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I'.'l 


O'.'e 


46.2 


0.8 


66.2 


0.8 


-0.087 
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0.009 
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1"0 


o'.'e 


46.6 


1.2 


61.9 


1.3 
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0.014 


0.014 


0.017 


0.049 


0.016 
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0.017 
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Major_2 


1'.'7 


o'.'e 


43.5 


1.1 


66.5 


1.1 


-0.081 


0.011 


0.016 


0.013 


0.035 


0.013 


-0.012 


0.013 


204 


Major_6 


1'.'7 


o'.'a 


44.3 


1.1 


64.6 


1.1 


-0.103 


0.011 


0.026 


0.014 


0.056 


0.013 


-0.020 


0.014 


204 


Major_5 


2'.'5 


C.'s 


46.3 


0.6 


61.3 


0.7 


-0.091 


0.007 


0.021 


0.009 
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Table 4 — Continued 











a" 






Ah3^ 


h4^ 


A/14'' 




Ahs'^ 




A/ie'^ 


S/N' Mask Name) 


32" 5 


2'.' 5 


-4.5 


1.9 


47.2 


2.4 


-0.059 


0.059 


-0.059 


0.071 


0.063 


0.065 


0.060 


0.071 


64 


40'.' 


5'.'0 


2.6 


2.9 


42.3 


3.9 


0.022 


0.117 


0.022 


0.142 


-0.020 


0.126 


-0.002 


0.141 


37 


50'.' 


5'.'0 


2.1 


2.7 


41.9 


3.6 


-0.040 


0.107 


-0.040 


0.131 


-0.004 


0.114 


0.062 


0.128 


42 


60'.' 


5'.'0 


-2.5 


4.0 


49.6 


5.3 


















30 


75'.' 


lO'.'O 


-2.2 


3.9 


39.0 


4.8 


















28 



Note. — The orientation of the axes has been defined as follows: 
major-axis runs NNW to SSE (0 — 160°), minor-axis runs ENE to 
WSW (0 = 250°). The systemie velocity of M32 is measured as 
— 199.7 lb 0.5km s~^ along the major-axis and —200.1 ± 1.8km s^^ 
along the minor-axis. 

^ Projected distance along the axis in areseeonds 
^ Spatial width of the bin in areseeonds 
^ Measured LOS velocity in km s~^ 

^ Scaled Poisson error in measured LOS velocity in km s~^ (where 
the scale factor / — 3.8 and 2.1 for the major- and minor-axis, respec- 
tively) 

^ Measured velocity dispersion in km 
f 

Scaled Poisson error in measured velocity dispersion in 
km s~^(where the scale factor / — 2.8 and 1.9 for the major- and 
minor-axis, respectively) 

^ Gauss-Hermite moments for spectra with S /N ^ 40 
^ Scaled Poisson error in the Gauss-Hermite moments (where the scale 
factor / — 1.9 and 2.1 for the major- and minor-axis, respectively 
^ Average S/N per pixel of the spectrum 

^ Mask names for the spectra that were not coadded (i.e. spectra 
located at |rM32| 2 ) 

2.3. An Integrated View of M32's Kinematics 

The results of M32's major- and minor-axis LOSVD 
analyses are shown in Figures [10] and 111) respectively. 
Our measurements push out much further that what has 
been observed to date. The velocity and velocity disper- 
sion profiles are measured out to ~ 230" along the major- 
axis and ~ 160" along the minor-axis using Gaussian 
fits to the resolved stellar population, with integrated 
light measurements reaching out to 85" along the major- 
axis and 75" along the minor-axis. Because the inte- 
grated light velocity measurement at 85" on the NNW 
major-axis of M32 is likely affected by M31 light (pulling 
the velocity data point more negative), we exclude this 
data point from our analysis in §[3l Measurement of the 
Gauss-Hermite moment profiles using integrated light ex- 
tends out to 70" along the major-axis and 50" along the 
minor-axis. The Gaussian fits to the resolved stellar pop- 
ulation extend well beyond the major-axis photometric 
distortion radius of 150" (minor-axis distance of 130") 
where surface photometry of M32 shows a sharp upward 
break in the surface brightness profile and elongation and 
twisting of the elliptical isophotes, distortions that may 
be the result of tidal interaction with M31. Our pro- 
files do not show evidence of sharp kinematical gradients 
across the distortion region. 

While the emphasis of our study is on the largely 
uncharted regions of M32, it is instructive to check 
whether our measurements agree with previous measure- 
ments of the inner regions. A co mparison of our in- 
tegrate d light measurements with Ivan der Marel et al.l 
(|1994b l) observations of M32's core is shown in Fig- 
ures [in] and [11] (right) . The profiles show good agree- 



ment, with minor differences attributable to spatial res- 
olution, slit position and template matching. Previous 
studies using integrated light to measure the mass of 
M32's central black hole have resulted in numerous de- 
tailed veloc ity measu rements of the core (iTonrvl [TOSTI: 
Dressier fc Richstoip 119881: Ivan der Marel et all 11994b ; 



Verolme et al.i i200a Ivan den Bosch fc de Zeeuwl I2010D 



The most radially extensive of these studies measures 
a two-dimensional mean velocity, velocity dispersion, 
and Gauss-Hermite moments h3 and fe 4 out to 30" us- 
ing wide field SAURON observations (iCappellari et al.l 
120071 Ivan den Bosch fc de Zeeuwl [20Tol ) . However, due 
to the large instrumental dispersion of the instrument, 
measurements for ve locity d ispersio n, /t3 and ft.4 may 
not be very accurate (jvan de n Bosch fc de Zeeu w 2010). 
Studies with higher instrumental resolution measure a 
velocity dispersion profile out to 13" along the major- 
axis and 22" along the miiior-axis ([Dressier fc Richstond 
119881 Ivan der Marel et al.l [r994bD . and Gauss-Hermite 
moments out t o 8" along the major-axis a nd 22" along 
the minor-axis (Ivan der Marel erallll994b[ ). 

M32's profiles appear smooth and symmetric. Our 
best-fit major-axis velocity profile rises steeply to a max- 
imum observed rotational velocity of Wrnax = 46 km 
at r = I'.'O (4pc). Our best-fit minor-axis velocity pro- 
file is relatively flat, with the small cusp seen at the 
center of the profile resulting from a 0"25 mis-centering 
of the longslit during observation. The observed cen- 
trally rising velocity dispersion profile peaks at CTmax ~ 
90 km s^^. These measurements are consistent with pre- 
vious ground-based observations. Observations of the 
core using a narrower slit with higher spatial resolution 
result in highe r measurements for veloci t y and veloc- 
ity d ispersion ([van der Marel et"aL[ Il997t [Joseph et alj 
[2OOII) . 

The systemic velocity of M32 is measured as — 199.7 ± 
0.5 km s^^ along the major-axis and —200.1 ± 1.8km s^^ 
along the minor-axis. These measurements are made by 
assuming point reflection symmetry of the velocity proflle 
at r < 10" along the major-axis, and outside the central 
cuspy region along the minor-axis. These measurements 
are consis tent with the Zwi cky Catalog value of —200 ± 
6 km s-i ([Falco et al.|[T999l) . 

3. DYNAMICAL MODELS 

3.1. Modeling Approach 

To interpret our new kinematical data we flt it with 
axisymmetric dynamical equilibrium models. The mod- 
els are constructed using Schwarzschild's orbit super- 
position technique. Over the past decade such mod- 
els have become a standard in the field, and they 
have been applied and tested repeatedly with different 
softw are im plemen tations on dat a for the galax y M32 
(e.g. Ivan der Marel et al.l Il998bl: iCretton et all 119991: 
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Figure 10. Best fit parameters to M32's major-axis (0 = 160° ) LOSVD from measurements of integrated ligfit (solid circles) and Gaussian 
fits to tlie resolved stellar population (open squares). Rows from top to bottom: M32's major-axis mean velocity, velocity dispersion, and 
tlie Gauss-Hermite moments /13, ^4, /15, and /ig profiles. Left column: Tfie full radial extent of eacli profile with la (68% confidence limit) 
error bars. Tlie integrated light measurement at — 85'.'0 (circled in red) is possibly affected by M31 contamination, thereby pulling the 
line-of-sight velocity point more negative. Right column: Comparison with van der Marel et al. (1994b:) measurements (open circles) shows 
good agreement. Error bars are not shown, but are generally smaller than the point size. The distance scale, assuming a distance to M32 
of 785 kpc, is 1" = 4pc. 

from our own only inside the central few arcsec, due to 
differences in spatial resolution. This is important for un- 
derstanding the black hole mass, but that is not a focus 
of the present study. At intermediate radii, where exist- 
ing datasets have some overlap with our own, the data 
are generally mutually consistent (see, e.g.. Figure ITU)) . 
Rather than to combine data from different authors, we 
therefore model here only the new Keck data, which span 
all observationally accessible radii along two axes. The 
quality of our fits discussed below, and their consistency 
with previous published work, lead us to believe that 
addition of integral fie ld data (available for < 28" from 
ICappellari et al.|[2006[ ) would not significantly alter the 
main results. 

For the modeling we use the software developed and de- 
scribed bv lvan der Marel et al.l (|1998b( ). with similar nu- 
merical parameter and grid settings as described therein. 
The models start with a luminosity density distribution. 



tended the modeling approach to include triaxial config- 
urations. However, they found that M32 is actually best 
fit by models that are close to axisymmetric and edge-on. 
We therefore restrict our analysis to such models here. 

Datasets of high spatial resolution obtained with 
the Hubble Space Telescope (|van der Marel erallll998al: 
iJoseph et al.ll200"l[) have demonstrated that M32 hosts a 
massive black hole in its center. Its mass is Mbh = (2.4± 
1.0) X 10^ Mq ()van den Bosch fc de Zeeuwll20lol . consis- 
tent with determinations by previous authors). The ob- 
servational setup for our new Keck data was optimized 
for spatial extent, and not spatial resolution. Rather 
than treat Mbh as a free parameter in our analysis, we 
therefore keep it fixed at (2.4 ± 1.0) x 10^ Mq in all our 
models. 

The data obtained by other authors generally differ 
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Figure 11. Same as Figure [Tol but for M32's minor- axis {(p = 250°). Th e slight mismatch seen at the center in the mean velocity profile 
(an apparent peak in our data and an apparent dip in Ivan der Marel et aL (1994b)) , and the fe3 and h5 profiles is an artifact that results 
from a slight mis-centering of the longslit in our observations and that of[yan der M arel et all l ll994bl ). 



with is transformed to a mass density under the assump- 
tion of a constant mass-to-light ratio M/ L. The luminos- 
ity density is chosen to fit the observed major axis sur- 
face brightness profile. For this we used the same V- band 
brightness profile as in Ivan der Marel eFall (|1998bh . but 
we tr ansformed this to the /-band using th e known V — I 
color (iLauer et al.lll998l: iTonr v et al."200I'). Slightly dif- 
ferent froml^^^^MareOLaL (1998b.) , we adopt a fore- 
ground extinct ion An = 0.35 and project e d axia l ratio 
q = 0.76 from Ivan den Bosch fc de Zeeuwl (|2010f) . This 
makes our M/L values directly comparable to theirs, 
since those authors also adopt the same distance as we 
do here. 

Given this approach, the M/ L is the only free param- 
eter to optimize the fit to the data. We construct mod- 
els with different M/L. The best-fit model is identified 
as the one that yields residuals with the overall mini- 
mum x^. We restrict our calculations here to models 
with a constant mass-to-light ratio for several reasons: 
(1) such models are easiest to calculate; (2) such mod- 
els have been found to adequately fit existing datasets 



for M32; (3) such models are a useful reference before 
considering more detailed modeling; and (4) such mod- 
els can be used in many cases to prove the presence of 
a dark halo (namely, if no constant mass-to-light ratio 
model can fit the data). 

3.2. Data-Model Comparison: Integrated- Light 

The predictions of the best-fit model are compared to 
the long-slit integrated light measurements in Figures [T^ 
and ll3[ for the major and minor axis respectively. 

The data-model comparison shows that integrated 
light kinematics for M32 are well-fit by an axisymmetric 
constant mass-to-light ratio model (with a central black 
hole). While the fit is not perfect (x^/Nuf = 2.35), all 
the trends in the data as function of radius are reason- 
ably well matched by the model. 

The finding that a constant mass-to-light ratio model 
adequately fits integrated-light kinematical data is con- 
sistent with what has been found by previous authors. 
However, our work extends this result to a radius that 
is three or more times larger than the region assessed 
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Figure 12. Data-model comparison for the long-slit data along the major axis. Data points with error bars are shown in black. Model 
curves shown in red display the predictions of the best-fit axisymmetric edge-on model. A constant mass-to-light ratio model with a central 
black hole provides an adequate fit out the the largest radii accessible with integrated light. 



by prior studies. This is a non-trivial finding, since one 
might have expected to start seeing the tell-tale signs of 
a possible dark halo at ~ Srog. But no such signs are 
readily evident. 

The (/-band) mass-to-light ratio of our b est-fit model 
is M/L = 1.24. Ivan den Bosch fc de Zeeuwl (pO lO) found 
acceptable (triaxial) models in the range M/L = 1.4 ± 
0.2. So while our best-fit value is lower than preferred by 
those authors, it is within the allowed range. The fact 
that, if anything, our fits over a much larger radial range 
yield a lower M/L is important. If in reality the velocities 
of M32 stars were elevated at large radii because of the 
presence of a dark halo, then fitting a constant mass-to- 
light ratio model should yield higher M/L values when 
data at increasing radii are included. We find instead 
the opposite. This strengthens the conclusion that the 
integrated-light data provide no strong indication for a 
dark halo in M32. 



3.3. Data-Model Comparison: Discrete Velocities 

The predictions of the best-fit model are compared to 
the V and a inferred from the discrete velocity measure- 
ments (blue open dots) in Figures [14] and [151 for the 
major and minor axis respectively. The long-slit data 
(black solid dots) are also included for comparison. 

A three- integral orbit superposition model has the free- 
dom to change its dynamical structure with radius. It 
is therefore worthwhile to consider for comparison sim- 
pler models that do not have this freedom. The green 
curves in Figure [T4| show the large radii predictions of 
two-integral models with a distribution function of the 
form f{E^ Lz). The Jeans equations of hydrostatic equi- 
librium can be explicitly solved for such models, making 
them a simple starting point for many analyses. They 
have been successfully applied to the case of M32 in many 
prior studies, and here w e have u sed the same software 
as in Ivan der Marel et all (jl994a[) to calculate their pre- 
dictions. The dotted green curves are for the same ge- 
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Figure 13. Same as Figure [121 but for the minor axis. The small central peak in mean velocity V is reproduced by a 0'.'25 perpendicular 
offset of the slit from the galaxy center. 



ometry and mass-to- light ratio as in our best-fitting orbit 
superposition model. We adopted a param eter k = 0.57 
(defined in e.g. Ivan der Marel et al]|l994al ) to split the 
azimuthal motion into ordered and random components. 

On the major axis, in the area around the ends of our 
long slits (40"-100" from the center), the two- integral 
model predictions are very similar to those from our best- 
fitting three- integral model (solid red curves). On the 
minor axis, the three-integral model predictions for the 
dispersion are higher than for the two-integral model. 
This may be because: (a) the best-fitting three-integral 
model has a dynamical structure that differs from a two- 
integral model; or (b) the three-integral model predic- 
tions are obtained from Gaussian fits to model LOSVDs, 
whereas the two-integral model predictions are true dis- 
persions. The long-dashed green curves in Figure [T5] 
show the minor-axis predictions for a two-integral model 
with a 21% higher value of M/L. These appear more 
similar to the three-integral model predictions. 

The purpose of the green two-integral curves in Fig- 



ures [14] and [15] is to show how the kinematics fall with 
radius in a constant mass-to-light ratio model in which 
the dynamical structure itself does not vary with radius. 
The gradient at large radii should not depend much on 
the dynamical structure itself, as long as it is indepen- 
dent of radius. The discrete velocity data points (blue 
points) at the large radii do not follow this nominal be- 
havior. In particular, the velocity dispersions in the four 
outermost data points are higher than expected. The 
two-integral models (using the higher M/L on the minor 
axis) predict on average cr — 24.7km s~^ at these radii. 
By contrast, the observed weighted average dispersion is 
33.0±3.5 km s~^, which is higher by a statistically signif- 
icant 2.4(7. This suggests an increasing M/L with radius 
(i.e., the presence of a dark halo), unless the dynamical 
structure of M32 changes with radius beyond the edges 
of our long slits. 

The three-integral orbit superposition modeling ap- 
proach automatically adjusts the dynamical structure as 
necessary to best fit all available data. The red dashed 
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Figure 14. Data-model comparison for tlie Gaussian parameters V and tr, along the major axis. Blue open data points are from fits 
to the discrete velocity measurement histograms. Solid black data points are from integrated-light measurements. Model curves shown 
in red display the predictions of the best-fit three-integral axisymmetric edge-on model. Model predictions for the integrated light data 
(corresponding to thin apertures) are connected by a solid red line. Model predictions for the discrete velocity data (corresponding to 
broad wedges, two on each side of the galaxy) are connected by a dashed red line. For comparison, dotted green curves at large radii show 
the major-axis predictions (not integrated over wedges) of a two-integral model with the same geometry and mass-to-light ratio. 



curves in Figures [U and [15] show the predictions thus ob- 
tained. For each side of the galaxy (major or minor axis, 
positive or negative radius) there are two data points. 
The predictions for these data points are connected by a 
straight line. These predictions correspond to averages 
over broad wedges on the sky, and not small apertures 
as was the case for the integrated-light predictions. This 
affects primarily the rotation velocity, which is smaller 
when averaged over a wedge than on the major axis it- 
self. 

Given the significant scatter between the discrete ve- 
locity data points, the three-integral model fits the data 
reasonably well. In particular, the model predictions sig- 
nificantly exceed the nominal model fall-off indicated by 
the two-integral models. For the four outermost data 
points the predicted dispersion is 32.7km s~^, consis- 
tent with the weighted average 33.0 ± 3.5 km s^^ of the 



observed values. For the four inner data points the pre- 
dicted dispersion does not fit the average of the observed 
values, 38.6km s"""^ versus 24.9±2.3km s~^, respectively. 
Apparently, a constant mass-to-light ratio model cannot 
simultaneously reproduce the low dispersions observed in 
the inner wedges of Figure [5l while also reproducing the 
higher dispersions observed in the outer wedges. 

The reason that the three-integral models predict 
higher dispersions than the two-integral models is due 
to a change in its dynamical structure. Inspection of 
the dynamical structure of the best-fitting three-integral 
model shows that it has increasing tangential anisotropy 
towards larger radii. This causes more motion to be ob- 
served along the line of sight direction. Moreover, mod- 
els with tangential anisotropy tend to have flat-topped 
LOSVDs. For such LOSVDs, the dispersion of the best- 
fitting Gaussian (which is the observed quantity) exceeds 
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Figure 15. Same as Figure [Ml but for the minor axis. The dotted green curves at large radii show the predictions of a two-integral 
model with the same geometry and mass-to-light ratio as the best-fit three-integral model (red). The long-dashed green curves show the 
predictions of a two-integral model with a 21% larger mass-to-light ratio. 



the true dispersion (|van der Marel k, Franxl[T993l ) . Both 
of these effects help the model to fit the observed disper- 
sions at large radii. 

To assess whether the tangential anisotropy of the best- 
fitting three-integral model is consistent with the data it 
is nec essary to measure th e shape of the LOSVD at large 
radii dCaroUo et al.l[T995l) . On average, the outer wedges 
in Figure [5] each have only some 28 observed M32 stars 
in them, much too little to reliable determine the Gauss- 
Hermite moments of the LOSVD. Nonetheless, some im- 
portant LOSVD shape information can be obtained from 
the data. Figure [16] shows the observed grand-total ve- 
locity histogram for the outer four wedges, with the M31 
contribution (held fixed at the values in Figure 6) sub- 
tracted. The red curve is the prediction of the best-fitting 
three-integral model. The curves have similar widths, as 
was already clear from the preceding discussion. How- 
ever, there are some subtle differences in the observed 
and predicted LOSVD shape. 



The vertical dashed curves in Figure [16] show the es- 
cape velocity of the model at i? = 200" in the equatorial 
plane (72km s"^). The wings of the predicted LOSVD 
fall to zero around this velocity. To achieve a signifi- 
cant dispersion, the model creats a relatively fiat-topped 
LOSVD within the regime bounded by the escape veloc- 
ity. By contrast, the observed LOSVD histogram has a 
narrower core, and broader wings. In particular, there 
are ~ 6 stars on the positive velocity side of the LOSVD 
that move faster than the model escape velocity. If these 
stars are bound to M32, then M32 must have a dark halo. 

However, alternative interpretations do exist. First, 
our simple model for the LOSVD of the M31 halo is a 
smooth Gaussian. If in reality the M31 halo has sig- 
nificant substructure (in velocity space) over the region 
surrounding M32, then this result may cease to be sig- 
nificant. It is then possible that the excess stars at w « 
95km s~^ in the M32 frame (i.e., Whci ~ —105km s^^) 
may simply be a co-moving group of M31 stars. Sec- 
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Figure 16. Observed grand-total velocity histogram for the outer 
four wedges in Figure with the M31 contribution (held fixed at 
the values in Figure[6]l subtracted. The red curve is the prediction 
of the best-fitting three-integral model. The horizontal axis shows 
the velocity in the M32 frame (djkji -|- 200 km s"-"^). Vertical dashed 
curves at ±72 km s~^indicate the escape velocity of the model at 
-R = 200" in the equatorial plane. The slight mismatch between 
the observed and predicted histograms may be interpreted either 
as evidence for an M32 dark halo, evidence for tidal perturbations 
in the M32 outskirts, or evidence for velocity substructure in the 
M31 halo. 

ond, it is possible that we have reached a regime in M32 
where tidal perturbations are playing a role. In this case, 
it would not be appropriate to interpret the excess stars 
in the context of an equilibrium model. The fact that the 
observed histogram in Figure [TB] is not symmetric (there 
are excess stars only on the positive velocity side) seems 
more consistent with either of these interpretations than 
with evidence for a dark halo. 

Other than the high-velocity tail of stars in Fig- 
ure 1161 our data-model comparison provides very little 
evidence that the kinematics of M32 might be affected 
by tidal perturbations or non-equilibrium dynamics. The 
integrated-light measurements inside 100" have a smooth 
behavior, and follow closely the predictions of our equi- 
librium dynamical models. The discrete measurements 
at larger radii show significant scatter, but this is likely 
due to measurement errors, and not tidally induced. The 
velocity dispersion increases from the inner to the outer 
wedges on Figure [5l but this increase is similar on all four 
sides of the galaxy, and therefore not easily attributed to 
tidal perturbations. 

In summary, by obtaining kinematics out to ^ Srcff 
we appear to have reached for the first time in M32 a 
regime where the observed velocity dispersion flattens 
out as a function of radius. This can be plausible inter- 
preted as evidence for a dark halo (although we have not 
actually constructed models with a dark halo explictly). 
On the other hand, a three-integral constant- Af/L model 
with increasing tangential anisotropy towards large radii 
can still fit all the data, apart from half-a-dozen stars at 
Vhc\ ~ —105 km s'-^. While these stars cannot be ex- 
plained by any equilibrium model without a dark halo. 



they may alternatively be due to substructure in the M31 
halo, or tidal perturbations in the M32 outskirts. 

We have not explicitly explored dynamical models that 
include a dark halo component. Given the added free- 
dom of a mass-to-light ratio that varies with radius, such 
models should certainly be able to fit the data. However, 
given the limited number of data points at the largest 
radii, it is unlikely that such models would be able to 
place strong constraints on the properties of the dark 
halo. 

4. SUMMARY & DISCUSSION 

We have presented mean velocity, velocity dispersion 
and higher order Gauss-Hermite moment profiles along 
M32's major and minor axes based on Keck/DEIMOS 
spectroscopic observations of the integrated light and 
the resolved stellar population. This study is the first 
to continuously transition between integrated light and 
the resolved stellar population in M32, or indeed in any 
galaxy. The kinematical profiles provide the most radi- 
ally extensive spectroscopic coverage for any cE galaxy, 
with measurements of the resolved stellar population ex- 
tending to a projected distance of ~ 8 r'j^ . 

We have constructed axisymmetric three-integral dy- 
namical models for M32 to interpret the new data. The 
integrated-light data out to ^ 100" show falling disper- 
sions with radius, which are well-fit by a constant M /L 
model. The discrete velocity data between 100-200" re- 
veal a regime where the observed velocity dispersion fiat- 
tens out as a function of radius. This can be plausible 
interpreted as evidence for a dark halo. However, a con- 
stant M/ L model can fit all the available data out to the 
largest radius probed, provided that there is increasing 
tangential anisotropy with radius in M32. The number 
of observed M32 stars at large radii is too small to di- 
rectly constrain the anisotropy. A small number of fast 
moving stars at large radii suggests that a dark halo may 
better explain the data than tangential anisotropy, but 
these stars may also be due to substructure in the M31 
halo, or tidal perturbations in the M32 outskirts. 

It has long been known that M32's isophotes undergo 
a sharp twist and increase in ellipticity coincident with a 
break in the surface brightness profile at 5 r'j^ , plan 



sibly as a result of tidal interaction with M31 (e.g. iKent 
1987tlGhoi et al.l 120021: [Johnston. Choi fc Guhathakurta 
20021 ). Our kinematical study, however, shows no cor- 
responding tidal signature across this region. On the 
contrary, M32's kinematics appear to be symmetric and 
regular out to the limit of our survey. The lack of a 
strong gradient in the velocity and velocity dispersion 
profiles across the radius at which isophotal distortion 
occurs does not necessarily rule out the tidal distortion 
hypothesis. An alternative explanation is that the ob- 
served isop hotal distortions i n M3 2 are intrinsic to its 
structure. iFasano fc Bmo^ (fT989t ) find an increase m 
ellipticity and isophotal twisting in the outer regions of 
about half of all isolated elliptical galaxies, and conclude 
that tides are unlikely to be the cause of the distortion 
seen in these galaxies. 

M32 is a compact elliptical galaxy, and such galaxies 
are normally found as satellites of more massive parent 
galaxies. As discussed in §[1] tidal stripping is a lead- 
ing theory for their formation. Our findings provide lit- 
tle support for a scenario in which tidal forces having 
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significantly altered the M32 structure at large radii for 
two different reasons, although our conclusions are not 
strong enough to rule the tidal stripping model out all 
together. First, the observed kinematics provide little 
direct evidence of being impacted by tidal effects. This 
appears surprising if the small size of M32 were directly 
attributable to such effects. Second, the discrete velocity 
measurements imply a relatively flat velocity dispersion 
profile at large radii. This suggests that M32 may well 
have a dark halo at radii of 200" and beyond, where there 
is very little luminous matter. Dark and luminous matter 
are both collision-less and subject to the same gravita- 
tional potential, and tidal stripping removes material at 
large radii first. So if a dark halo was left prominent at 
radii of 200" and beyond, then it may be difficult to ex- 
plain why the luminous body of the galaxy got stripped 
down to a radius Tcff ~ 30". 

Neither the observed kinematics nor the dynamical 
modeling implications about the presence of a dark halo 
provide evidence that tidal forces have significantly al- 
tered the M32 structure at large radii. These findings 
alone do not rule out tidal models for cE formation but 
instead suggest that tidal forces may not be as important 
in shaping cE galaxies as often argued. 
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